NIFA Visualizations

library(haven)
library(ggplot2)
library(plotly)

## 
## Attaching package: 'plotly'

## The following object is masked from 'package:ggplot2':
## 
##     last_plot

## The following object is masked from 'package:stats':
## 
##     filter

## The following object is masked from 'package:graphics':
## 
##     layout

categ <- read_dta(file = "ByCateg_v5.dta")

nutrients <- read_dta(file = "foodcode_nutrient_density.dta") 

categfull <- categ %>%
  filter(!is.na(price)) %>%
  filter(categ_name != "Uncat")

Initial Data Visualizations

ggplot(data = categfull, mapping = aes(x = price, y = densprot)) +
geom_point(aes(color = broadcateg_name, size = tspend), 
           alpha = 0.5) +
  labs(title = "Protein on a Budget?", x = "Price/1000 kcal", y = "Protein/1000 kcal") +
  scale_size(guide = 'none') +
  scale_color_discrete(name = "Food Category")

ggplot(data = categfull, mapping = aes(x = price, y = denssodi)) +
geom_point(aes(color = broadcateg_name, size = tspend), 
           alpha = 0.5) +
  labs(title = "Sodium on a Budget?", x = "Price/1000 kcal", y = "Protein/1000 kcal") +
  scale_size(guide = 'none') +
  scale_color_discrete(name = "Food Category")

ggplot(data = categfull, mapping = aes(x = price, y = dens_add_sugars)) +
geom_point(aes(color = broadcateg_name, size = tspend), 
           alpha = 0.5) +
  labs(title = "Sugar on a Budget?", x = "Price/1000 kcal", y = "Protein/1000 kcal") +
  scale_size(guide = 'none') +
  scale_color_discrete(name = "Food Category")

ggplot(data = categfull, mapping = aes(x = price, y = densfibe)) +
geom_point(aes(color = broadcateg_name, size = tspend), 
           alpha = 0.5) +
  labs(title = "Fiber on a Budget?", x = "Price/1000 kcal", y = "Protein/1000 kcal") +
  scale_size(guide = 'none') +
  scale_color_discrete(name = "Food Category")

First Plotly Iteration

Protein

# Legend
leg <- list(
  bordercolor = "#E2E2E2",
  borderwidth = 0.5)

# Protein 
prot <- plot_ly(
  categfull, x = ~price, y = ~densprot,
  # Hover text:
  text = ~paste("<b>Type:</b>", categ_name, "<br><b>Description:</b>", categ_des),
  color = ~broadcateg_name, size = ~tspend, opacity = 0.8
)

prot %>%
  layout(title = "Protein",
                     xaxis = list(title = "Price/1000kcal"), yaxis = list(title = "Protein(gm)/1000kcal"),
                     margin = list(l = 123), legend = leg) 

Fiber

# Legend
leg1 <- list(
  bordercolor = "#E2E2E2",
  borderwidth = 0.5)

# Protein 
fibe <- plot_ly(
  categfull, x = ~price, y = ~densfibe,
  # Hover text:
  text = ~paste("<b>Type:</b>", categ_name, "<br><b>Description:</b>", categ_des),
  color = ~broadcateg_name, size = ~tspend, opacity = 0.8
)

fibe %>%
  layout(title = "Fiber",
                     xaxis = list(title = "Price/1000kcal"), yaxis = list(title = "Fiber(gm)/1000kcal"),
                     margin = list(l = 123), legend = leg1) 

Plotly Visualization with Dropdown Feature

#Plot with Dropdown
prot2 <- plot_ly(
  data = categfull, x = ~price, y = ~densprot,
          color =  ~broadcateg_name, size = ~tspend,
          alpha = 0.7,
          width = 900,
          height = 500,
          hoverinfo = "text",
          text = ~paste0("<b>", categ_name, "</b><br>", "<b><i>Description: </b></i>", categ_des, "<br><b><i>Daily Spending/Capita: </b></i>", round(tspend, 3)),
    #ifelse(categfull$categ_name == categfull$categ_des, (~paste0("<b>", categ_name, "</b><br>", "<b><i>Description: </b></i>", categ_des, "<br><b><i>Total Spending/Capita: </b></i>", tspend)), (~paste0("<b>", categ_name, "</b><br>", "<br><b><i>Total Daily Spending: </b></i>", tspend))),         
          type = "scatter",
          mode = "markers"
          ) %>%
  layout(
    title = list(text = "<b>Finding Affordable Nutrient Sources</b>", x = .75),
    xaxis = list(title = "Price($)/1000kcal"),
    yaxis = list(title = ""),
    margin = list(b=150), ##right margin in pixels
    annotations = 
           list(x = 1.32, y = .4, #position of text adjust as needed
                text = "*Datapoint size reflects \nfood category's total \ndaily spending", 
                showarrow = F, xref='paper', yref='paper', 
                xanchor='right', yanchor='auto', xshift=0, yshift=0,
                font=list(size=12, color="grey")), 
    legend=list(title=list(text='<b> Broad Food Category </b>')),
       updatemenus = list(
        list(
          y = .55,
          x = -.05,
        type = "list",
        label = 'Category',
      buttons = list(
         list(method = "restyle",
              ##list(method = "update",
                  ##args = list("y", list(~densprot)),
          args = list(list(y = list(categfull$densprot))),
               #args = list("y", list(categfull$densprot)),
               label = "Protein(gm)/1000kcal"),
          list(
           ## method = "update",
           ##  args = list("y", list(~densfibe)),
            method = "restyle",
               args = list("y", list(categfull$densfibe)),
               label = "Fiber(gm)/1000kcal"),
          list( 
            ## method = "update",
            ## args = list("y", list(~dens_add_sugars)),
      method = "restyle",
               args = list("y", list(categfull$dens_add_sugars)),
               label = "Added Sugar(gm)/1000kcal"),
          list(
            ## method = "update",
            ## args = list("y", list(~denssodi)),
            method = "restyle",
              args = list("y", list(categfull$denssodi)),
               label = "Sodium(mg)/1000kcal")))))

prot2